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' The electron energy distribution index, p, is a fundamental parameter of the process by 

i which electrons at shock fronts are accelerated to rclativistic speeds and by which they 

radiate, via synchrotron emission. This acceleration process is applicable to a myriad of 
astronomical jet sources such as AGN, X-ray binaries and gamma-ray bursts (GRBs), 
as well as to particle acceleration in the solar wind and supcrnovae, and the acceleration 
^ ' of cosmic rays. The accurate measurement of the distribution of p is of fundamental 

importance to differentiate between the possible theories of electron acceleration at any 
rclativistic shock front; there is division as to whether p has a universal value or whether 
it has a distribution, and if so, what that distribution is. 

Here one such source of synchrotron emission is examined: the blast waves of GRB 
, afterglows observed by the Swift satellite. Within the framework of the GRB blast wave 

f ^ ■ model, the constraints placed on the distribution of p by the observed X-ray spectral and 

temporal indices are examined and the distribution parametrized. The results show that 
the observed distribution of spectral indices is inconsistent with an underlying distribu- 
tion of p composed of a single discrete value but consistent with a Gaussian distribution 
centred at p = 2.4 and having a width of 0.6. This finding disagrees with theoretical 
work that argues for a single, universal value of p, but also demonstrates that the width 
|TN| ' of the distribution is not as wide as has been suggested by some authors. 

Keywords: Gamma-ray burst: general; Radiation mechanisms: non-thermal; Acceleration 
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1. Introduction 

The afterglow emission of Gamma- Ray Bursts (GRBs) is generally accepted to orig- 
inate from external shocks when a collimated ultra- rela tiyis tic jet ploughs into the 
circumburst medium, driving a shock front ahead of iiUJ2l. The shock accelerates 
charged particles in the jet or ambient medium to the relativistic speeds required so 
that they may radiate electromagnetically via synchrotron emission. This accelera- 
tion mechanism - common to many astronomical jet sources such as AGN, X-ray 
binaries and GRBs (as well as particle acceleration in the solar wind and super- 
novae, and the acceleration of cosmic rays) - is thought to be Fermi diffusive shock 



1 



5, 2011 0:25 WSPC/INSTRUCTION FILE CurranT 



2 P. A. Curran 

acceleration^ due to the passage of an external shocli^^ after which the energy of 
the electrons, E, follows a power-law distribution, N(E)dE oc E~ p dE, with a cut- 
off at low energies. Fundamental to this is the electron energy distribution index, p, 
the characteristic parameter of the process, which is dependent only on the under- 
lying micro-physics of the acceleration process. Some (semi-)analytical calcul ations 
and simulations indicate that there is a nearly universal value of ~ 2.2 — 
though other studies suggest that there is a large range of possible values for p of 
1.5-M. Observationally, different methods have been applied to samples of BATSE, 
BeppoSAX and Swift GRBs which reached the conclusion t hat the observe d range 
of p values is not consistent with a single central value of jJ 10 | H | 12 | 13 | 14 | an j ^ ne 
latter three showed that the width of the parent distribution is a,, ~ 0.3 — 0.5. In 

I -i r I 

these proceedings I summarize our previous worW^ in which we interpret a much 
larger and, statistically, more significant sample of Swift observed GRB afterglows, 
to constrain the electron energy distribution, p. I further compare the distributions 
found from spectral and temporal indices, as well as those found for other sources. 



2. Method and Results 

The method, detailed in Ref. 1151 is to constrain the electron energy distribution 
index, p, from the values of the X-ray spectral indices, /3x or /3, observed by the 
Swift XRT up to July 200£P^. p is derived from the spectral index as opposed to the 
temporal index because for a given spectral index, assuming the asymptotic limit, 
there are only two possible values of p depending on whether the cooling frequency, 
v c , is less than or greater than the X-ray regime, vx, while for a given temporal 
index there are m ultip le possible values which are model dependent (e.g., the simple 
blast wave modet^, or modifications thereof). In accordance with synchrotron 
emission predicted by the blast wave model, the behaviour of the unabsorbed X-ray 
spectrum is described as a single power law, where the flux goes as F v (y) oc . 
Under the standard assumptions of slow cooling and adiabatic expansion of the blast 
wave, the electron energy distribution index is related, in the asymptotic limit, to 
the spectral index by either p — 2/3 (y c < vx) or p — 2/3 + 1 {v c > vx) (e.g., 
Ref. [17j) : implying a difference between the slopes of the two spectral regimes of 
A/3 = 0.5. The regime probability, X, is defined as the probability that the cooling 
frequency is less than the frequency of the X-ray regime (i.e., v c < vx) and 1 — X 
is the probability that v c > vx (Figure [T]). The case where the cooling frequency 
may be passing through the X-ray regime is neglected (since there is no sign of 
spectral evolution in the sample) as are the cases where the peak frequency, v m , 
or self absorption frequency, z^ a , are greater than the X-ray regime as these are not 
observed in late time afterglows. To parameterise the underlying distribution of the 
electron energy distribution index, p, a maximum likelihood fit is used to return 
the most likely parameters of the assumed underlying model, the errors on which 
are estimated via a Monte Carlo error analysis. Another Monte Carlo analysis tests 
the probability that the observed distribution of spectral indices could be obtained 
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Fig. 1. The left curve shows the probability distribution function of the underlying value of p, 
while the right, double-peaked curve shows the transformation of that into X-ray spectral index 
space. The distance between the peaks is equal to the difference between the spectral indices either 
side of a cooling break (A/3 = 0.5), while the relative heights of the two peaks are dependent on 
the relative probability that the X-ray regime falls above or below the cooling frequency. 

from an underlying distribution of p described by the most likely parameters. 

There are two underlying models, or hypotheses, regarding the data that we 
wanted to test: that the observed distribution of spectral indices, ft, can be obtained 
from an underlying distribution of p composed of i) a single discrete value and 
ii) a Gaussian distribution. We founcP^ that the observed distribution of spectral 
indices are inconsistent with the first hypothesis but consistent with the second, 
a Gaussian distribution centred at p = 2.36 and having a width of 0.59 (Figure [2] 
Left). Furthermore, if one accepts that the underlying distribution is a Gaussian, the 
majority (> 94%) of GRB afterglows in our sample have a cooling break frequency 
less than the X-ray frequency. 

3. Discussion 

To further test the c onc lusion regarding the distribution of p, the values of p pub- 
lished for 191 GRBa^ are tested in a similar way to above using Monte Carlo 
analyses. It is found that these data (Figure [2] Right) are well fit by a Gaussian 
distribution of p, centered at p — 2.36 and with a width of a — 0.36. It is impor- 
tant to note that this is a totally independent data set from that described above - 
based on temporal fits of XRT light curves as opposed to spectral fits - and as such 
offers an independent test of our results for a similar Swift-XRT GRB sample. The 
excellent agreement between these two independent data sets offers confirmation of 
the results for the distribution of p, as well as our conclusion that the majority of 
GRB afterglows in the sample have a cooling break frequency less than the X-ray 
frequency, since both are interdependent in our study. 
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Fig. 2. Le/t: Normalised histogram of the spectral data (301 measurements) overlaid with high- 
resolution normalised histograms of the synthesized data sets (10 4 X 301 data points) from the 
most likely parameters of a single discrete p (narrow distribution) and a Gaussian distribution of 
p (wide distribution). Note that there is only one obvious peak on the spectral index distribution 
while from the transformation of figure [T] two were expected. Right: Normalised histogram of the 
191 values of p published in Ref. 1181 overlaid with the probability density function of p = 2.36 and 
u = 0.36, which is consistent with that data. 

Th is result confirms the results from previous small-sample GRB afterglow 
studies^HSISl as regards the non-universality of p, the central value at p ~ 2.0—2.5, 
and the width of the distribution of a p ~ 0.3 — 0.5. However, our results are based 
on a sample of bursts an order of magnitude larger than these studies and a statis- 
tical approach to the position of cooling frequency relative to the X-ray is taken, by 
using the regime probability, X. The result is also consistent with (semi- ) analytical 
calculations and simulations that indicate a value of p ~ 2.2 — 2.4SEEI Though it 
rules out the suggestion of a single, universal value, the result demonstrates that the 
distribution is not as wide as has been suggested by some authors (1.5 < p < 4^. 

Given that these results suggest the value of p is not universal, it may be pos- 
sible that p changes suddenly or evolves gradually with time or radius even in a 
single event as environmental shock parameters (e.g., magnetic field, ambient den- 
sity) change or evolve. It is also possible that different components of a structured 
jet, multi-component jet or jet-cocoon could have different values of p. If either are 
the case our derived values of p should be considered time averaged values of the 
parameter, though a change or evolution of p should be observable as a change or 
evolution of the synchrotron spectral index, ft, and no significant example of such 
an evolution has been observed in GRB afterglows. 

It is very clear from the data (Figure [2] Left) , as well as the fitting, that there 
is no secondary peak as had been expected (Figure Q]) and this implies that the 
majority of GRB afterglows in the sample (> 94%) have a cooling break frequency 
less than the X-ray frequency. While we found a distribution of p that concurs 
with previous work, this is somewhat anomalous, as previous multi-band (optical - 
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X-ray) studieJUHSHH have shown that the cooling break frequencies of a number 
of GRB afterglows are greater than the X-ray frequency. The upper limit derived 
implies that this is true in only ~ 20 of 301 GRBs in the sample and while this 
is a discrepancy from the ratio observed in the multi-band studies (5 out of 10^1 
and 2 out of it is not seriously so, given the extremely low number statistics of 
those studies. In fact, if a sample of 6 random bursts is created from the statistical 
distribution of 301, there is a non-negligible (~ 5%) chance that 2, or more, bursts 
would have a cooling frequency greater than the X-ray frequency, consistent with 
the aforementioned Swift based studjEl, which investigates the same population. 

The lack of a second peak in the data, if not intrinsic, could be due to a very 
wide underlying distribution of p which would blur out the double peak structure; 
however this requires a far wider distribution than is observed. It might also be 
explained by a bimodal distribution of p where the mode from which p is drawn 
is correlated with position of cooling break but this is a highly convoluted and 
unlikely scenario. Another possible explanation is that instead of p being related 
to the spectral index through one of the two relationships presented in section 
there is only a single relationship or a continuous range of relationships. This 
is certainly possible, most obviously caused by the fact that the canonical GRB 
spectra does no t have sharp breaks between the regimes but is, instead, smoothly 
transitioning^! This would cause a continuous range of relationships which would 
manifest as a single peaked spectral index distribution centered about the average 
of the two relationships, /3w(p- 0.5)/2, and which would, in the case of our data 
imply that p w 2.9 with a width of <r p ~ 0.6. However, importantly the cooling 
frequency evolves in time so any such spectra due to a smooth transition between 
regimes would likewise evolve and no such significant evolution has been observed 
in these GRB afterglow!^, leading us to conclude that this is unlikely to be an 
effect and that the single peak structure is, in fact, intrinsic. 



3.1. Comparison with other sources 

A similar, though smaller sample, effort has been made to derive the distribution 
of p in blazars and pulsar wind nebulae^l which found that p « 2.4 (a p » 0.4) 
and p « 1.7 (a p ~ 0.6) respectively. But the uncertainty of these authors' methods 
is demonstrated by the fact that it was also possible that the distribution could 
be described asp ~ 3.2 (a p » 0.2) and p w 2.0 (a p « 0.2) respectively, when 
the distribution was forced to be as narrow as possible. Though, no such attempt 
has been made to statistically analyze X-ray binaries, the observed radio spectra of 
these sources have spectral indices of order 0.5-0.7, corresponding to p ~ 2.0-2.3, 
while a small sample of FR I active galaxies have been treated to a non-statistical 
analysis^' which found a narrow range at 2.1. It is planned to expand the above 
described statistical method to these other sources to obtain accurate measurements 
of the distribution in the different classes. These may then be used as a touchstone 
for theories and simulations and will allow proposed models to be tested. 
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4. Conclusions 

These results for Swift observed GRBs show that the observed distribution of spec- 
tral indices are inconsistent with an underlying distribution of p composed of a 
single discrete value but consistent with a Gaussian distribution centred at p = 2.4 
and having a width of 0.6. This finding disagrees with the theoretical work that 
argues for a single, universal value of p, but also demonstrates that the width of the 
distribution is not as wide as has been suggested by some authors. Due to the lack 
of accurately measured distributions of p for other sources, an accurate comparison 
with those sources is not possible, however, it is planned to expand this statisti- 
cal method to other sources to obtain accurate measurements of the distribution 
in the different classes, which may then be used as a touchstone for theories and 
simulations. 
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